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Abstract 

Axial-vector interactions of neutrinos and axions with a medium are based on 
structure functions which cannot be calculated in simple perturbative terms for nu- 
clear densities. We use the SN 1987A neutrino signal duration to estimate the range 
of allowed neutrino opacities of a supernova (SN) core. We perform numerical sim- 
ulations of the evolution from the hot, lepton-rich post-collapse stage to the cold, 
deleptonized neutron star and compare the predicted neutrino signal characteristics 
with the SN 1987A data. Besides the neutrino opacities we vary the mass, temper- 
ature profile, and equation of state of the initial model. Within plausible limits of 
these quantities the neutrino opacity cannot be much smaller than its "standard" 
value. This limits the width of the spin-density structure function to much smaller 
values than indicated by a perturbative calculation, implying that bremsstrahlung 
processes such as axion emission probably saturate at around 10% nuclear density. 
A substantial suppression at yet larger densities that might have been expected does 
not seem to occur. 

PACS numbers: 14.60.Lm, 14.80.Mz, 95.30.Cq, 97.60.Bw 
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1 Introduction 



Aside from gravitational waves neutrinos provide the only known possibility to obtain 
direct information about the physical processes in matter at the extreme conditions of 
density and temperature that occur at the center of type-II supernova (SN) explosions, in 
colliding or coalescing binary neutron stars, or in collapsing white dwarfs. A type-II SN 
explosion marks the violent final stage in the life of a massive star with a main sequence 
mass 8M0 ^ M ^ 30 M©. It is associated with the disruption of the star and an outburst 
of light with a total energy of about erg. Most of the stellar mantle and envelope is 
ejected with a kinetic energy of the order of 10^^ erg, forming the ultimately diffuse SN 
remnant. The progenitor star's core with a mass of 1-2 M0 is thought to form a compact 
object, i.e. a neutron star or possibly a black hole. 

Core collapse events and merging binary neutron stars are powerful sources of neutri- 
nos which are far more efficient at carrying away energy from the dense inner regions than 
photons for which the opacity is about 20 orders of magnitude larger. 99% of the gravita- 
tional binding energy (a few 10^^ erg) of the neutron star forming in a SN are emitted in 
neutrinos and antineutrinos of all flavors. They carry away lepton number and energy, thus 
driving the evolution from the hot, lepton-rich collapsed stellar core to the cold, neutron- 
ized remnant. In addition, recent progress in our understanding of SN physics indicates 
that neutrinos play an active and crucial role for different processes during type-II SN 
explosions. Neutrino energy deposition in the layers between the newly formed neutron 
star and the SN shock front is thought to yield the energy for ejecting the stellar mantle 
and envelope [1]. This energy deposition creates a "hot bubble" outside of the protoneu- 
tron star which has comparatively low density, but quite high temperatures [2]. Also, the 
transfer of energy from neutrinos to the material near the surface of the protoneutron star 
causes a continuous flow of matter into the hot-bubble region. This neutrino-driven wind 
may well be the long-sought site for the r-process nucleosynthesis [3]. 

The measured SN 1987A neutrino signal at the Kamiokande II and 1MB water Che- 
renkov detectors [4] crudely confirms theoretical expectations. Notably, the duration of 
the signal of several seconds precludes efficient cooling via "invisible channels" such as the 
emission of right-handed neutrinos or axions. Therefore, the SN 1987A signal has been 
used to set a large number of constraints on a variety of novel particle physics models 
and phenomena [5]. Using SN 1987A and perhaps a future galactic SN as a "particle 
physics laboratory" along these lines, as well as a proper understanding of the SN explosion 
mechanism and r-process nucleosynthesis all require a reliable calculation of the "neutrino 
light curve" from a core collapse SN. 

To this end the transport of neutrinos through the stellar medium has to be simulated 
in numerical models of SN explosions and of the protoneutron star evohition, ideally by 
solving the Boltzmann transport equation. In practice, an approximate treatment based 
on the corresponding moment equations is usually employed. A necessary input are the 
dynamical structure functions of the stellar medium which account for the interactions 
with the neutrinos. For instance, the neutral-current scattering ^{ki) — > i'{k2) is com- 
pletely described by a function W{ki, k2) which gives the probability for a transition from 
a neutrino with four momentum ki to one with k2. For most practical applications this 
function was calculated on the basis of a simple cross section for the scattering process 
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^{ki) + N ^ N + v{k2) on free nucleons A^. In the nonrelativistic limit nucleon recoils can 
be neglected; the scattering cross section is then found to be cr = ((jp/47r) (Cy + 3C|^) £\ 
with the Fermi constant Gp, the vector and axial- vector coupling constants^ Cy and Ca, 
and the initial neutrino energy ei. The scattering cross section and the neutrino opacities 
are thus found to be dominated by axial-vector (spin-dependent) interactions. The same 
holds true for cliarged-current reactions of the form + n -H- p + e~ and I^e + p -H- n + e"*" . 

In a recent paper Raffelt and Seckel [6] have questioned the validity of this approach for 
the conditions of a SN core. Simply put, they argued that the collisions among nucleons 
which interact by a spin-dependent force cause the nucleon spins to fluctuate fast relative 
to other relevant time scales. For neutrino scattering such a time scale is the inverse of 
the energy transfer Ae = 62 — £i- Neutrinos cannot "resolve" processes which happen at a 
faster rate and so, if the spin of a given nucleon flips several times during the period (Ae)"^ 
the neutrino will see an approximately vanishing nucleon spin. For the vector current, a 
similar effect does not exist because in the nonrelativistic limit there are only small nucleon 
velocity fluctuations. This difference between the vector and axial-vector interactions is 
also illustrated by the absence of neutrino pair bremsstrahlung N + N^N + N + V + u 
from the vector current, while there is a large contribution from the axial vector. In this 
case neutrino pairs are emitted by the fluctuating nucleon spins which act as a time- varying 
source for the neutrino axial- vector current. 

In an extreme scenario, the axial-vector cross section could be almost entirely sup- 
pressed at high densities, leaving only the vector-current interaction as an opacity source. 
At first sight this picture appears ruled out by the good agreement between the SN 1987A 
neutrino observations with standard predictions. In a preliminary study [7] it was found, 
however, that this agreement might be deceiving in the sense that a modified neutrino spec- 
trum and a modified emission period could conspire such as to leave the observed event 
numbers and signal durations relatively unchanged. Presently we formalize this study in 
order to understand the impact of modified (reduced) neutrino opacities as predicted by 
Raffelt and Seckel [6] . We perform numerical simulations of the neutrino cooling of nascent 
neutron stars for modified neutrino opacities and compute the expected neutrino signals in 
the Kamiokande II and 1MB detectors. Asking for consistency between these signals and 
the observed ones we attempt to infer limits on the possible amount by which the effective 
neutrino cross sections could be reduced in the stellar plasma. 

The paper is organized as follows. In Sect. 2 we summarize the considerations in favor 
of a possible suppression of the axial- vector weak currents at SN conditions. We motivate a 
simple formula for the suppression effect as a function of density and temperature. A para- 
metric ansatz allows us to switch continuously between full suppression and no suppression 
of axial-vector currents in our numerical simulations of the neutrino-driven cooling of pro- 
toneutron stars (Sect. 3.1). Aside from modified neutrino opacities for neutral-current 
and/or charged-current processes, we consider models with different neutron star masses, 
with different versions of the equation of state (EOS) of matter around and above nuclear 

^In vacuum, the values are C{? = -1, = 1 - AsiiT^ @w = 0.07, Q « -1.15, and « 1.37. The 
previously used values = — = 1-26 had to be modified because some of the nucleon spin appears 
to be carried by strange quarks. (For references see [6].) In a bulk nuclear medium these coefficients are 
likely suppressed, in analogy to the suppression of the charged current where Ca is suppressed from its 
vacuum value of 1.26 to about 1.0 in a nuclear medium. 



3 



density, and with different initial temperature profiles. Section 3.2 contains a comparison 
of expected and observed neutrino signals, and Sect. 4 closes the paper with a discussion 
and interpretation of our results. 



2 Neutrino scattering rate 

2.1 The dynamical spin-structure function 

We are mostly concerned with the neutral-current axial-vector scattering rate of neutrinos 
on the nuclear medium. As a simple model we consider a medium consisting of only 
one species of nucleons with a number density equivalent to the baryon density n^. The 
differential scattering cross section per nucleon for the transition i'{ei) ^'{£2} can then 
be expressed in the form 

d£^ ~ ' ^ ' 

where Gp is the Fermi constant and Ca the neutral-current axial-vector coupling constant 
with the values for protons and neutrons discussed in the introduction. 

In the limit of nonrelativistic nucleons, and ignoring nucleon recoils as well as collective 
spin motions, one finds that S'o.(Ae) is the spin-density dynamical structure function in 
the long- wavelength limit (e.g. [6], [8]) 



SJuj) = lim / die'"*(s(i,fc) •s(0,-fc)). (2) 

SrzB fe^-o J-^ ^ ' 

Here, Si = tp^cTitp (i = 1, 2, 3, Pauli matrix cr^, nucleon Pauli spinor tp) is the spin-density 
operator. A calculation of this correlator (and more complicated ones in a mixed medium of 
several types of constituents) is the main problem at treating neutrino transport correctly. 

For a medium of nucleons which do not interact with each other, i.e., in the limit of a 
very dilute medium one easily finds 

S^{lo) = 2TrS{uj). (3) 

This form refiects that neutrinos do not change their energy in collisions with isolated 
"heavy" nucleons (no recoil!). The total axial- vector scattering cross section is then cta = 
SClGlej/in. 

In the long-wavelength limit the structure function governs all processes which rely 
on a coupling to the nucleon spin such as the axial-vector vV emission or absorption or 
the emission or absorption of axions. Axions are hypothetical pseudoscalar particles which 
couple to nucleons by virtue of the axial-vector interaction Lagrangian (C/2/„) V-'T/iTs'*/' 
where C is a model-dependent dimensionless number of order unity and /„ is the axion 
decay constant or Peccei-Quinn scale, a quantity with the dimension of an energy. If 
axions can escape freely from a neutron star so that one may ignore re-absorption and 
Bose-stimulation effects, the energy-loss rate per unit volume is given by 
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In the dilute limit where Sfj{u)) 2nS{Lo) the axion emission rate vanishes. 

When nucleon-nucleon interactions are not ignored, the lowest-order contribution to 
axion emission is bremsstrahlung N + N N + N + a. An explicit calculation was 
performed by Brinkmann and Turner [9] in what corresponds to the long- wavelength limit. 
They modelled the nucleon-nucleon interaction by a one-pion exchange potential with 
massless pions. Their result can be brought into the form 

S^{uj) = r^uj-^s{uj/T) (5) 

where T is the temperature and represents something like a nucleon spin rate of change 
due to collisions. It is found to be [6] 

= T - Ty^mf ~ '%o I T ) 

where a^^ = (/2mjv/m^)^/47r ~ 17 with / ^ 1.05 is the "pion-nucleon fine structure 
constant" and po = 3x 10^'' g cm~"^ is the nuclear density. The factor a;~^ is easily shown to 
represent the nonrelativistic propagator of the intermediate nucleon in the bremsstrahlung 
graph [6] . It can also be shown to follow from a classical treatment of the bremsstrahlung 
process [10]. 

The remaining dimensionless function s{x) of the dimensionless energy transfer x = 
Lo/T represents the nucleon phase space volume in A^A^^-collisions. For positive energy 
transfers (energy given to the neutrino in a collision) it is found to be 

/■oo 

six)=e-- / dye-yi\x\y + y'Y/'. (7) 
Jo 

It has the asymptotic behavior s(0) = 1 and s{x ^ 1) = e^*(a;7r/4)^/^. For negative energy 
transfers one finds that the factor e~^ must be replaced by 1. Then 

S^iuj) = S^i-uj)e--/^ (8) 

as required by detailed balance. 

The total cross section based on Eq. (1) together with Eq. (5) diverges because of the 
uj~^ behavior of the structure function. Raffelt and Seckel [6] proposed to modify this 
divergence by a Lorentzian cutoff uj~^ {u>^ + F^/4)~^ with a suitable energy scale F. Its 
magnitude is dictated by the normalization requirement 

/ + 00 
duj S^{uj) = 2^ (9) 
-oo 

which can be shown to follow from the definition of if there are no correlations between 
the spins of different nucleons in the medium. For <^T one finds F = F^-, i.e., the struc- 
ture function S(j{uj) is essentially a Lorentzian of width F^r which for small Fq^ approaches 
the (5-function already found for a gas of non-interacting nucleons. 

The normalization condition illuminates the intuitive fact that the total neutrino scat- 
tering rate on nucleons is not increased just because the nucleons interact with each other 
before and after the interaction with the neutrino. However, these "upstream" or "down- 
stream" collisions do allow for the transfer of energy between nucleons and neutrinos that 
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would not be possible otherwise because of the smallness of nucleoli recoils. It is also 
clear that one cannot conceptually separate the elastic single-nucleon scattering process 
v + N ^ N + u from the inelastic process u + N + N — > N + N + v because any nu- 
cleon with which a neutrino can interact is known with certainty to have emerged from 
a previous iViV-coUision, and with certainty will interact again. Therefore, the emerging 
picture is that Su{uj) is a narrow S-\ike function in a dilute medium, and that it develops 
"wings" as the medium becomes denser. The "wings" are probably well determined by the 
lowest-order bremsstrahlung calculation while the region around a; = is dominated by 
multiple scattering effects. 



2.2 Sa-{ijj) at high density 

As long as ^ T the picture of neutrino scattering on quasi-free nucleons and the 
bremsstrahlung version of S(^{uj) can be reconciled with each other without much problems 
by the above Lorentzian ansatz for the low-o; behavior. This "dilute limit" corresponds 
physically to the picture that a nucleon spin flips only rarely relative to the time scale set 
by inverse thermal frequencies of the system. In a SN core, however, one is confronted 
with the opposite limit 3> T, i.e., nominally the nucleon spin flips many times within 
one oscillation period of, say, a thermal neutrino. 

In order to remind the reader of the typical physical conditions present in a young SN 
core we show in Fig. 1 the proflles of various relevant parameters for our reference case, 
model S2BH_0 (see Sect. 3), about 1 s after the numerical simulation of the neutrino cooling 
was started. This stage of the evolution corresponds to the situation in a protoneutron 
star 1-2 s after its formation. The degeneracy parameter of neutrons is typically between 
2 and 4, i.e., neutrons are mildly degenerate. However, for most purposes a value r/ ?a 3 
marks the borderline between degenerate and nondegenerate conditions whence we think 
it is justifled to use the nondegenerate limit for our exploratory study in this paper. The 
parameter T^/T is found to be typically of order 30, i.e., r„ ^ T. 

It is not at all obvious what Sa-{io) will look like for such conditions. Besides the detailed- 
balance condition Eq. (8) which must apply on fundamental grounds one may assume that 
the normalization condition Eq. (9) remains valid (absence of spin-spin correlations) , and 
one may assume that for w ^ F^- the bremsstrahlung result Eq. (6) remains valid. Together, 
these three conditions imply that Sfj{Lo) is a broad and shallow function as opposed to a 5- 
function in the dilute limit. As a generic example for the shape of Scr{uj) in the high-density 
limit we require the asymptotic high-o; behavior Sa^{u!) = T^^ (7r/4)^''^7o- x^'^/^ e^* where 
X = Lo/T. Moreover, we speculate that for a; ^ F^- the structure function is essentially a 
constant. Then, an ansatz for its shape is 



Here, Scr{x) is the dimensionless function TScr{Tx). For 7o- ^ 1 the normalization integral 
is dominated by negative energy transfers whence (3 = {7r/27)'^^'^y^. This example illus- 
trates the implications of the normalization condition in conjunction with the perturbative 
calculations which flx the high-a; asymptotic form. 
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FIG. 1: Physical parameters of our reference model S2BH_0 after about Is of neutrino cooling 
and deleptonization, corresponding to the situation in a protoneutron star about 1-2 s after 
stellar core collapse. The plots in the left column show profiles of density p (top), temperature T 
(middle), and the number fractions of electron neutrinos, Y^, electrons, Y^, and leptons (electrons 
plus electron neutrinos), Yl (bottom) as functions of the enclosed baryonic mass Mb. The right 
column of figures displays the Fermi momenta of neutrons, protons, electrons, and neutrinos 
(top), the corresponding degeneracy parameters of n and p as ratios of the Fermi energies to the 
temperature (middle), and the parameter = F^/T of Eq. (6) (bottom). 
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The impact of this high-density behavior of S„{u)) on the neutral-current neutrino 
opacity is crudely estimated by the inverse mean free path, averaged over a thermal en- 
ergy spectrum of the initial neutrino. Moreover, all expressions become much simpler if 
one replaces the Fermi-Dirac occupation numbers with the Maxwell-Boltzmann expression 
g-£i/T. |.]^g resulting error is small for nondegenerate neutrinos. The relevant quantity is 
then 

(An \ 1 
-j) = ^l ^^^l dx^xlxle-'^^S^ix^-x,), (11) 

where Xq^ = {SCj^GpUB / 4:Tr) 24 T^. One energy integration can be done explicitly, leaving 
one with an integral over the energy transfer alone {x = lo/T) 

1 f°° 

F = — dx{\2 + Qx + x'^)S„{x). (12) 

IZTT Jo 

This expression is unity for elastic scattering, i.e., when Sa{x) = 2'n-S{x) or at least a 
narrow, normalized function. 

For 7^ » 1 one finds with Eq. (10) F (135\/3/27r2) 7^ ^ 12 7^ ^ Therefore, the 
total neutrino scattering rate drops precipitously for 7o- larger than a few. In order to 
interpolate to the low-density regime where F = 1 one may use F = (1 + 7^/12)"^ where 
7o- = Tf^/T was given in Eq. (6). With ~ 20 in the SN core, the axial- vector neutral- 
current scattering rate would be entirely suppressed. 

While this radical conclusion follows from our three assumptions concerning S^^, one or 
more of them may not apply. We are most suspicious of the assumption that the width of 
Sfj{Lo) is indeed given by To- of Eq. (6) as this quantity essentially represents the spin rate 
of change due to nucleon-nucleon collisions. As colliding nucleons themselves are subject 
to other collisions, their spins fluctuate as well, leading to an averaging effect which may 
well reduce below its perturbative value. It does not seem plausible that the spins will 
fluctuate much faster than thermal frequencies given by the energy levels of the system. 
Thus, one may speculate that does not grow much broader once its width is of order 
T. This sort of behavior would imply, for example, that the axion emission rate reaches 
some maximum level at which it saturates without dropping at yet higher densities. The 
axial-vector neutrino scattering rate would be suppressed, but perhaps not entirely. 

The main goal of our present study is to find out if the neutrino signal of SN 1987A 
allows us to distinguish between such speculations, and notably, if it allows us to exclude 
the most radical scenario of a complete suppression of the axial-vector scattering rate. 

Keeping in line with these exploratory goals of our study we implement modified neu- 
trino opacities by the prescription C\ FC\ with 

^ = (1 - ^) + if^ (13) 

where 

. 1 (^A"" ( P V^ lOMeV ^ 

'=Y2[t) "i 3xl0i3gcm b J (^^j " ^^'^ 

Here, a = 1 corresponds to the full suppression effect while a = to no suppression at all. 
For completeness we also consider a < 0, implying enhanced cross sections. 
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While the discussion of Raffelt and Seckel [6] focussed on neutral-current processes one 
may speculate that similar reduction effects would pertain to charged-current processes. 
Therefore, we will also investigate the case where the same axial-vector reduction factor 
applies to both types of reactions. 

3 Neutrino cooling of protoneutron stars 

During the first seconds of its life the newly formed neutron star evolves from a hot, lepton- 
rich state to a cold, neutronized object, a process driven by the loss of lepton number and 
energy via the emission of several 10^^ neutrinos and antineutrinos of all flavors [11, 12]. 
In the present work we are interested in simple characteristics of the neutrino "light curve" 
during this quasistatic "Kelvin-Helmholtz cooling phase" such as the total energy emitted 
in neutrinos, their average energies, and the expected number of events and their time 
distribution at the Kamiokande II and 1MB detectors. In particular, we focus on the 
variation of these characteristics when the axial-vector neutral-current scattering rates of 
neutrinos are assumed to be suppressed to a certain degree. As a matter of comparison we 
also vary other aspects of uncertain input physics such as the equation of state (EOS) and 
the baryonic mass and temperature profile of the initial protoneutron star model. 

3.1 Numerical simulations of protoneutron star evolution 

The set of initial models used as an input into our evolution simulations was obtained by 
scaling parameters such as mass and temperature of a protoneutron star model from a 
core collapse calculation [13] at a stage about half a second after the formation of the SN 
shock. The EOS of matter around and above nuclear density po is a finite-temperature 
extension of a zero-temperature EOS developed by Glendenning [14, 15], which optionally 
allows us to choose a gas of neutrons, protons, electrons, and muons ("EOS A") or to 
include additional baryonic states like hyperons and A-resonances at densities above ^ 2po 
("EOS B"). For more details about this EOS and the structure of Wilson's post-collapse 
model, for information about the numerical procedure and the implementation of EOS and 
neutrino physics, and for a description of the typical evolution of a cooling and neutronizing 
protoneutron star, see Ref. [12]. 

The nomenclature of our models was adapted from the paper of Keil and Janka [12]. 
For example, the letters "SBH" ("Small-EOS B-Hot") mean a protoneutron star model 
with a mass smaller than the maximum, stable neutron star mass for EOS B and with an 
initial central temperature which yields a gravitational mass only little below the baryonic 
mass, thus ensuring that the model reproduces the marginally bound state of a stellar core 
shortly after collapse. Correspondingly, "SAH" denotes a model for EOS A. Systematic 
extensions to these names indicate variations of model properties made in the present work. 
The sample of investigated models can be grouped in the following way. 

1. Models evolved with different neutrino opacities: The degree of the suppression of 
axial-vector neutral and charged currents is measured by the free parameter a in 
Eq. (13). The chosen values a = x.yz with 0.00 < x.yz < 1.00 are indicated by 
the extensions ''jryz" in the model names, a = 0, i.e. the standard case of the 
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neutrino cross sections, is simply marked by "_0", negative values — 1 < a < 
(corresponding to enhanced cross sections) are marked by the suffix "m" behind 
the numbers. The specification 'Lxijzt" means models where the critical density in 
Eq. (14) was changed from 3 x 10^'^g/cm'^ to 10^''g/cm\ When only the neutral 
currents for neutrino-nucleon interactions are modified, we add an "n" at the end of 
the names. 

2. Models with different masses, initial temperatures, and EOS's: From low to high 
baryonic mass the models are labelled sequentially as "SOB...", "SIB...", "S2B...", 
etc. The models of the sequence S#BH_0 were constructed such that the ratio of 
the initial gravitating mass Mg,o to the baryonic mass Mb is constant. Models with 
a higher initial temperature, i.e. bigger initial gravitational mass, have an asterisk, 
"*" , at the end of their names, even higher temperatures are marked by "**" . When 
EOS A instead of EOS B was used, "B" in the model specification is replaced by an 
"A" . 

In Table 1 we present an overview over our sample of models. We list the characteristic 
parameter combination of each model, i.e. the baryonic mass Mb, the initial gravitational 
mass Mg,o, the initial baryon number density n^o at the center of the star, the initial 
central temperature Tcfl, and the value of a. In addition, the symbols used in the figures 
are indicated. 

Results of the simulation runs are also collected in Table 1. The deleptonization time 
scale ^di measures the time until the protoneutron star has reached the deleptonized state, 
in which the electron lepton number adopts its lowest possible value, corresponding to a 
chemical potential //^^ = for electron neutrinos, is the duration of the cooling process 
by neutrino emission and is defined as the time it takes until the temperatures have dropped 
below ^ 3 MeV in the entire star. E^tot denotes the total energy (redshifted for an observer 
at infinity) that is radiated away during time tc\, while tgo gives the time interval after 
which 90% of this energy have left the star. The energy emitted in Vg, which may be 
absorbed on protons in the water Cherenkov detectors of the 1MB and Kamiokande II 
experiments and which are responsible for most of the signals there, is listed in the column 
^^EpJ\ We characterize the major differences in the overall behavior and the evolution 
of different models by these simple parameters. For a description of the typical cooling 
history of a nascent neutron star we refer the reader to the paper by Keil and Janka [12]. 

3.1.1 Modified neutron star parameters 

From the data listed in Table 1 and plotted in Fig. 2 we learn that the time scales t^i, 
tc\, and tgo as well as the energies Etot and E^,^ increase with growing neutron star mass. 
The increase of the time scales is moderate and essentially linear for masses 1.3 M© ^ 
M ^ 1.6 M0, but becomes dramatic for larger masses. The cooling time scale and the 
time vary by more than a factor of 3; ^ci changes from about 28 s for a small star with 
Mb = 1.30 M0 to nearly 100 s for the 1.75 M0 model, which is very close to the mass 
limit of stable neutron stars for EOS B, M^^^ = 1.77 M0. The occurrence of hyperons 
in a large part of the star has a strong infiuence on the cooling evohition [12]. A more 
massive protoneutron star loses more gravitational binding energy, the total energy release 
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FIG. 2: Cooling time scale t^, deleptonization time scale t^i, time interval igo of 90% of the energy 
loss by neutrino emission, and time intervals <kii and tiuB of 90% of the predicted neutrino counts 
in the Kamiokande II and 1MB detectors, respectively, as functions of the baryonic mass of the 
protoneutron star. Different plot symbols correspond to the different models listed in Table 1. 
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FIG. 3: Predicted numbers of neutrino counts TVkh and A^imb in the Kamiokande II and 1MB 
detectors, respectively, and the corresponding average energies {£[,^)kii and (ep^)iMB of captured 
electron antineutrinos vs. the baryonic mass of the protoneutron star. Different plot symbols 
correspond to the different models listed in Table 1. 
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being roughly proportional to the increase in baryonic mass with a rate of change of about 
0.4 X 10^^ erg per 0.1 M0. This dependence becomes less steep towards the high-mass end 
of our model sample where the production of hyperons consumes a considerable amount 
of energy and general-relativistic redshift effects gain more importance. 

The effects of redshift and of hyperon production at the high densities in massive pro- 
toneutron stars can clearly be seen from a comparison of model S4BH_0 with S4AH_0 
which have the same baryonic and the same initial gravitational mass. In model S4AH_0, 
however, hyperons are absent (EOS A); it releases a larger amount of energy because 
general-relativistic effects are stronger by about 30% in model S4BH_0, where hyperoniza- 
tion affects as much as the inner 1.2 Mq. Also, deleptonization and cooling proceed sig- 
nificantly more slowly in model S4BH_0 because the occurrence of hyperons in the central 
region leads to a more compact core with higher densities and higher temperatures and 
correspondingly higher neutrino opacities. The peak temperatures in model S4BH_0 reach 
65 MeV, in model S4AH_0 they climb up to only 55 MeV. The cold, deleptonized star 
S4BH_0 has a central baryon density of 1.037 fm^'^ and a radius of 10.7 km, whereas the 
final density of S4AH_0 is = 0.519 fm~^ and the radius 12.2 km. 

Higher initial temperatures and thus, larger initial gravitating masses (models S2BH_0*, 
S2BH_0**, and S3BH_0*) lead to an additional energy release of roughly 0.23 x lO'^'^erg per 
0.01 Mq. In all cases, the energy radiated away as VeS is about 0.160-0.162 (slightly less 
than 1/6) of the total energy loss. 

3.1.2 Modified neutrino opacities 

Model S2BH_0 with a baryonic mass of 1.53 Mq is used as a reference star for varying the 
neutrino opacities. Changing the "suppression parameter" a (Sect. 2.2) from —1 to 1, we 
find a monotonic decrease of the time scales td, t^i, and tgo. The cooling and deleptonization 
of the protoneutron star become very fast as a — > 1. shrinks from a value of 108 s for 
a = —1 (model S2BH_100m) to 41s for the case of standard opacities, a = 0, and drops 
by another factor of 9 to only 4.5 s for nearly complete suppression of axial-vector weak 
currents in model S2BH_095 with a = 0.95 (Table 1, also Figs. 4 and 5). Interestingly, the 
relative amount of energy leaving the star as electron antineutrinos decreases monotonically 
when the stellar matter is assumed to be more transparent to neutrinos. 

If the charged-current interactions are left unchanged, the modification of the neutral- 
current reactions has a somewhat lesser impact than before, in particular the variation of 
the deleptonization time scale is damped. The overall picture, however, remains unchanged. 
Replacing the "critical density" in the denominator of Eq. (14) by the value lO^^g/cm^ 
has some influence, too, but again the general trends prevail. These results indicate that 
neutral-current neutrino-nucleon scattering reactions play the dominant role as an opacity 
source. This is understood because all types of neutrinos (/Vg, 77e, t'^, V^, and Vt) are 
affected, and because electron neutrino and antineutrino absorptions are strongly blocked 
due to fermion degeneracy. Moreover, most of the matter of the nascent neutron stars is at 
densities above 10^^ g/cm^, the surface layers yield only a minor contribution to the optical 
depth. 
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FIG. 4: Cooling time scale t^], deleptonization time scale idh time interval tgo of 90% of the energy 
loss by neutrino emission, and time intervals iKii and tiuB of 90% of the predicted neutrino counts 
in the Kamiokande II and 1MB detectors, respectively, as functions of the suppression parameter 
a of axial- vector neutrino-nucleon interactions. The plot symbols correspond to the models listed 
in Table 1. 




FIG. 5: Enlargement of the region 0.4 < a < 1 of Fig. 4. 
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FIG. 6: Predicted numbers of neutrino counts N^u and TVimb in the Kamiokande II and 1MB 
detectors, respectively, and the corresponding average energies (£i;e)Kii and {ep^)iMB of captured 
electron antineutrinos vs. the suppression parameter a of axial-vector neutrino-nucleon interac- 
tions. The plot symbols correspond to the models listed in Table 1. 

3.2 Comparison with SN 1987A neutrino signal 

In order to judge the degree of parameter variations which is compatible with the SN 1987A 
neutrino signal, we computed the expected neutrino signals of our models, in particular the 
signal durations ixii and iiMB, the expected event numbers A^kii and A^imb, and the average 
energies (ei^Jkii and (£,5^)imb of electron antineutrinos absorbed in the Kamiokande II and 
1MB experiments, respectively. The corresponding data are listed in Table 1. A/^kii and 
NiMB represent the total numbers of counts predicted for each model, ^kii and ^imb give 
the time intervals in which 90% of them occur. The source was assumed at a distance of 
50 kpc, the approximate distance of the Large Magellanic Cloud. Details of the procedure 
to evaluate our models to obtain these data can be found in the appendix of the work by 
Keil and Janka [12]. 

3.2.1 Modified neutron star parameters 

The average energies of the neutrinos measured in both detectors reveal an amazing stabil- 
ity against differences of the masses of our models (see Table 1 and Fig. 3) . This indicates 
some form of self-regulation of the position of the neutrino sphere, which is the layer below 
the surface of the star where neutrinos decouple from the stellar background. Due to the 
higher core temperatures of more massive stars the neutrinos are more energetic. However, 
the neutrino opacity of the stellar matter also increases because the interaction rates of 
neutrinos depend on the particle energies and thus on the temperature. As a result the po- 
sition of the neutrino sphere moves further out into the cooler regions closer to the surface 
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of the star. 

The number of events and signal duration are increasing functions of the neutron star 
mass (Table 1, Figs. 2 and 3). Deviations from the essentially linear dependence become 
visible again for our most massive model S4BH_0 with Mb = 1.75 M©. The signal duration 
t-Kii reflects the steep growth of the deleptonization and cooling times for the heavy stars 
(Fig. 2). In contrast, ^imb tends to decrease (Fig. 2), and iVxn and A^imb are nearly inert 
to the larger energy loss of model S4BH_0 compared with model S3BH_0 (Mb = 1.65 M0), 
see Fig. 3 and Table 1. This behavior is caused by the strong gravitational redshift of the 
massive model S4BH_0 which is very compact because hyperonization softens the EOS in 
a large part of the star. In particular, the 1MB experiment with its high energy threshold 
(about 20 MeV) is sensitive to the corresponding suppression of the high energies in the 
neutrino spectrum. A comparison between models S4BH_0 and S4AH_0 confirms this 
interpretation and again demonstrates the considerable effects that result from disregarding 
hyperons in the much stiffer EOS A. 

Our initial models were constructed "by hand" rather than obtained from a self- 
consistent collapse calculation. This entails an uncertainty which may be estimated from 
the changes of the calculated detector responses for models with modified initial tem- 
peratures (gravitating masses), i.e. models S2BH_0, S2BH_0*, and S2BH_0**, and models 
S3BH_0 and S3BH_0*. The signal durations are completely unaffected while the number 
of counts scatters by 2-4 units in the investigated cases. 

Provided "standard" neutrino opacities and EOS B yield a reasonable description of 
the physics at the extreme conditions in SN cores, our sample of cooling simulations with 
modified neutron star masses allows for the conclusion that the SN 1987A neutrino pulse 
most likely indicates a protoneutron star with a baryonic mass between 1.5 and 1.6 Mq. 
We caution, however, that in the present study we did not intend to find the best-fitting 
model for the neutrino source in SN 1987A. 

3.2.2 Modified neutrino opacities 

Reduced or enhanced neutrino opacities cause clear changes of the neutrino signature in the 
Kamiokande II and 1MB experiments. The predicted signal durations follow the trends of 
td\, tci, and tQQ. For the very opaque case (a = —1, model S2BH_100m) Kamiokande II and 
1MB measure events over periods of 37 s and 19 s, respectively (Table 1 and Fig. 4), while 
for a nearly complete suppression of the axial-vector contributions to neutrino-nucleon 
interactions (a = 0.95, model S2BH_095) the signal durations drop to about 2 s in both 
detectors (Fig. 5). 

Interestingly, the expected number of counts and also the mean energies of neutrinos 
captured in the detectors show the inverse behavior. When the value of a is changed 
from a = — 1 to 0.95, iVxn increases from 9 to 21, and A^imb from 3 to 26. This strong 
sensitivity can be understood by the spectral characteristics of the emitted neutrinos. 
Neutrinos decouple from the stellar medium at an optical depth of about unity. For a low 
neutrino opacity of the neutron star matter they leave the star from deeper inside, i.e. from 
layers of higher densities and temperatures. Therefore, these neutrinos are more energetic 
as is evident from the monotonic rise of (e^jKn and (epJiMB with higher values of the 
suppression factor a (Fig. 6). 
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Modifying only neutral-current interactions while keeping the charged currents fixed 
(models marked with an "n" at the end of their names) leads to changes of the predicted 
neutrino signals relative to our standard case S2BH_0 which are only slightly less dramatic. 
For model S2BH_090r with a higher value of the critical density, we find, in contrast, an 
even stronger increase of the mean energies of the captured neutrinos and, correspondingly, 
of the number of detector counts than for model S2BH_090 or model S2BH_090n. This is 
caused by the stronger heating of the outer regions of the star where the cooling of the 
low-opacity inner part is slowed down by the less transparent surface layer. 

The opposite trends of the different signal parameters (decreasing cooling and delep- 
tonization time scales, increasing count rates and mean neutrino energies) allow us to 
exclude a certain range of a as being highly unlikely. The extremely short signal durations 
and the very large numbers of energetic neutrino events predicted for models S2BH_090, 
S2BH_095, S2BH_090n, and S2BH_090r seem to rule out a complete or nearly complete 
suppression of axial- vector current interactions^. Neither statistical fluctuations of the 
data nor theoretical uncertainties like the EOS of supranuclear matter or the baryonic 
and gravitational mass of the protoneutron star seem to be able to account for the signal 
characteristics that obtain in the limit a — > 1. 

3.2.3 Modified neutrino spectra 

The total number of events observed in the 1MB detector is rather sensitive to the high- 
energy tail of the I7e-flux because of the relatively high energy threshold of 1MB. In the 

protoneutron star cooling calculations presented here, neutrino transport was treated by an 
equilibrium-difi'usion method [12], implying that the spectral I?e-distribution is essentially 
a Fermi-Dirac function with a degeneracy parameter r; = 0. A detailed Monte Carlo 
treatment of neutrino transport [16] reveals, however, that the I7e-spectrum is probably 
much better approximated by a Fermi-Dirac function with r] = 2-3. If the average Ve- 
energy is kept flxed this implies that both the low- and high-energy parts of the spectrum 
are depleted relative to intermediate energies, i.e., the spectral distribution is "pinched". 
This effect would cause a substantial depletion of the observable number of events at 1MB. 

In order to study the impact of this effect quantitatively we took the T^g-flux and average 
I7g-energies {cp^} of our reference case (model S2BH_0) and calculated the expected neutrino 
counts, the average e"^-energies in the Kamiokande II and 1MB detectors, and the expected 
signal durations as functions of rj. We stress, again, that (epj was held fixed at a given 
time, not the spectral temperature, which was adjusted for a given r] to yield the same 
{sp^). In Fig. 7 we show the relative change of Nku, -/Vimb, (ee+)Kii, {£e+)iMB, *kii, and 
^iMB- Only iViMB responds sensitively to r]. For the standard range t] = 2 — 3 one obtains 
a reduction between 30% and 50%. 

The uncertainty introduced by the unknown amount of pinching, the possibility that 
the neutrino spectra may have been modified by oscillations (e.g. Ref. [17]), and the small 

^The Kamiokande II signal consists of eight events within about 2 s and another three events between 
about 9 and 12 s after the first one. Our conclusions are based on taking the late-time events as part 
of the neutrino signal, i.e., the long time gap of about 7s at Kamiokande is interpreted as a statistical 
fluctuation. 
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FIG. 7: Impact of the "spectral pinching" on the total numbers of events in the Kamiokande II 
and 1MB detectors (A^kii, -^imb), on the average energies of the observed positrons from the 
reaction Vg + p ^ n + e'^ , and on the expected signal durations iKii and iiMB- 

integrated number of events in each detector together reveal that the total numbers of 
events, and especially Nimb, are rather poor measures to learn much about the properties 
of the neutrino signal. The most sensitive observables remain the durations of the signals 
registered at the two detectors. 

4 Discussion and conclusions 

Recently, Raffelt and Seckel [6] conjectured that the axial- vector part of weak neutral 
currents — possibly also of weak charged currents — in neutrino-nucleon interactions 
should be (partially) suppressed at the high densities present in SN cores. This effect would 
be a consequence of fast fluctuations of the nucleon spins due to frequent collisions. In the 
present work we investigated the implications of this possibility in the astrophysical scenario 
of neutron star formation. For that purpose we performed numerical simulations of the 
neutrino-driven cooling and deleptonization of newly formed neutron stars and determined 
the changes associated with a reduction of the neutrino cross sections. We computed the 
expected neutrino signals for the Kamiokande II and 1MB detectors. Consistency between 
these theoretical predictions and the SN 1987A neutrino observations can be used to obtain 
limits on the acceptable range of neutrino opacities. 

Our simulations showed that the duration of neutrino emission from the nascent neu- 
tron star and the corresponding signal durations at the detectors are drastically reduced 
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in the limit of a (nearly) complete suppression of the axial-vector current interactions. 
Interestingly, due to the increased transparency of neutron star matter for neutrinos this 
is accompanied by a strong rise of the number of neutrino counts and, in particular, of the 
average energies of the absorbed electron antineutrinos. These diverging trends lead to a 
clear discrepancy of the characteristics of the theoretical neutrino signals with the observed 
ones. 

Testing the sensitivity of the signal against variations of unknown properties of the 
nascent neutron star (mass of the protoneutron star, initial temperature profile, EOS at 
supranuclear densities) we found that the associated uncertainties cannot compensate the 
changes from extremely suppressed neutrino opacities. In particular, it seems impossible 
to construct a scenario where both of the diverging trends are simultaneously compen- 
sated. For example, a decrease of the signal time due to reduced neutrino opacities of the 
protoneutron star could be moderated by a higher mass of the star. However, in that case 
the number of events will rise even more. Therefore, our numerical results demonstrate 
that a strong suppression of the axial- vector part of neutrino-nucleon interactions can be 
firmly excluded. Being more quantitative, it appears very unlikely that suppression factors 
of a ^ 0.5 do occur. 

Concerning some hypothetical mechanism that might lead to an amplification of neu- 
trino-nucleon interactions at high densities we are not able to draw equally definite conclu- 
sions. Even though the signal durations become appreciably longer when the parameter 
a gets negative, this could be compensated, at least in part, by the faster cooling of a 
smaller protoneutron star. The mean energies of the captured neutrinos do not display 
a similar sensitivity for a < and thus cannot serve to put strong limits on an accept- 
able enhancement of neutrino opacities. Nevertheless, our results indicate that a slight 
decrease of (ejjJiMB is already sufficient to lead to a considerable decrease of the number 
of neutrino counts in the 1MB experiment. This adds to the already suppressed number 
of neutrino events from a less massive neutron star. Therefore, an enhancement of the 
neutrino opacities by more than 50% can be considered as quite improbable and by a fac- 
tor of 2 or more seems to be excluded. However, we stress that the count numbers in the 
detectors, in particular in the 1MB detector, are rather poor measures of the properties of 
the neutrino emitting star. On the one hand they are sensitive to the uncertain amount 
of spectral pinching (Sect. 3.2.3), on the other hand they turned out to change noticeably 
with the temperature of the star (Sect. 3.2.1) or by additional neutrino emission due to 
mass accretion onto the newly formed neutron star [18]. 

Therefore, the observed neutrino signal indicates that the neutrino opacities cannot be 
too different from what a naive cross section calculation fov v + N ^ N + v yields which, 
in turn, is based on a spin structure function of the form Sa{(^) = 2tvS{uj). This does 
not imply, however, that a naive perturbative calculation of all microscopic reaction rates 
involving the axial-vector current was essentially appropriate. A naive calculation of the 
bremsstrahlung process N + N N + N + V + u or N + N —?• N + N + a (axion a) entails 
that Sa{uj) should be a broad function of width ^ T where the "spin ffuctuation rate" 
To- was defined in Eq. (6). Our findings indicate that Sfj{Lo) is a function with a width not 
larger than a few T; this conclusion depends on the normalization condition of S^, i.e., the 
assumption of uncorrelated spin motions. Put another way, our findings appear to imply 
that the spin fiuctuation rate in a nuclear medium is of order T and not as large as given 
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by a naive perturbative estimate. 

Even though the neutrino opacity appears to be similar to what a calculation based on 
free nucleons yields and what was used in practical applications in the past, the axial- vector 
structure function is, no doubt, very broad rather than a (5-function. Therefore, in a single 
collision the energy of the neutrino typically changes by a large amount. This is in contrast 
to the standard picture where the neutral-current collision on "heavy" nucleons (small 
recoil) conserves the neutrino energy. According to Eq. (6) the width of Sa{uj) is around 
3T already at 10% nuclear density so that the picture of energy-conserving neutral-current 
collisions breaks down at relatively low densities, whether or not the structure function 
becomes even broader at higher densities. (Our study suggests that it does not become 
much broader.) We have not addressed this issue in our present work, even though it may 
have an impact on the predicted neutrino spectrum and notably on the "pinching effect" 
found in the Monte Carlo simulations of Janka and Hillebrandt [16]. 

The bremsstrahlung emission of neutrino pairs or axions is given by an integral of the 
type Eq. (4); it involves a factor e^'^/^ which is included in Sa-{io) for positive energy 
transfers. For a normalized Su{uj) this integral is maximal if Su{lo) is a broad function 
with a width of a few T. As our results seem to indicate that S^iuj) does not grow much 
broader than a few T one is led to conclude that the bremsstrahlung emission rate saturates 
at its maximum possible value which is reached at around 10% nuclear density, without 
dropping at yet larger densities. Granting this, the axion emission rate from a SN core is 
most realistically estimated by its value at the given temperatures and about 10% nuclear 
density. For typical axion models this means that nia ^ 10~^ eV is excluded, i.e., previous 
SN 1987A axion bounds [5] are not diminished as much as one would have feared on the 
basis of the extreme suppression scenario discussed in Sect. 2. 

While the SN 1987A neutrino signal indicates that the neutrino opacities are not too 
far from their "standard" values, their actual magnitudes have not been computed on the 
basis of first principles and so, they essentially remain adjustable parameters to account for 
a SN neutrino light curve. Because the opacities are dominated by the axial-vector neutral- 
current interactions with the nuclear medium one can estimate a plausible functional form 
of the spin-structure function S(^{uj) which, in turn, governs quantities such as the axion 
emission rate. Therefore, even though the axion emission rate can not be calculated by 
simple perturbative means one can estimate its magnitude from the "observed" neutrino 
opacities. As a result, a SN neutrino signal remains a useful tool to constrain particle 
physics phenomena based on the axial-vector coupling such as the axion-nucleon coupling. 
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'ABLE 1: Compilation of protoneutron star models with their characteristic parameters, the baryonic mass Mb, and the initial values of gravitational 
lass Mg,0 5 central baryon density n^o, and central temperature T<._o- In addition, the suppression factor a of axial- vector contributions to neutrino- nucleon 
iteractions and the plot symbol are given for each model. The neutrino-driven evolution of the protoneutron stars is described by the time scales of the 
eleptonization process, tdi, and of the cooling, igo is defined as the time interval after which 90% of the total energy emitted in neutrinos, -Btot, have left 
le star. £^tot is the energy measured by an observer at infinity. The fraction of this energy which is released as electron antineutrinos is given by Ej^^ . As for 
le predicted neutrino signals in the Kamiokande II and 1MB detectors, we list the integrated numbers of events in the detectors, Akh and Nimb, the time 
itervals iKii and iiMB in which 90% of these events occur, and the average energies of the captured electron antineutrinos, (e^jKii and (epJimb; respectively. 
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